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Abstract 



We compute the polarisability of the nucleon to leading order in chiral perturbation 
theory. The contributions from kaons and baryon resonances as intermediate states are 
included in addition to the contribution from pions and nucleons that had been previously 
computed. The isoscalar operators are dominated by the infrared behaviour of pion loops 
giving rise to a l/m^j^ coefficient. In contrast, the isovector operators are dominated by 
loops involving kaons, giving a 1/mfc coefficient, and further demonstrates that the strange 
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quark is an important component of the nucleon. In addition, the inclusion of the decuplet 
of baryon resonances as intermediate states substantially modifies the result found from 
the octet baryons alone for the isoscalar polarisability. 
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The electromagnetic polarisability of the nucleon is a unique probe of the structure 
of the nucleon. Like the charge radius and anapole moment, it is dominated by long- 
distance physics in the form of pion and kaon loops and diverges in the chiral limit as 1/m 
(m being the mass of the pion or kaon). Recent experimental measurements of the both 
the proton [|1| and neutron polarisabilities have generated increased theoretical effort 
to understand at a more fundemental level what these measurements tell us about the 
structure of the nucleon. Several theoretical approaches have been used to estimate the 
polarisability of the nucleon. These include; dispersion relations 0] constituent quark 
models 0-0, bag models [|TD[-[0, pion cloud models the Skyrme model []TB|-[]T^ 
and chiral perturbation theory (for a review see |]2l|). In the context of chiral 

perturbation theory, previous computations of the polarisability at one loop have included 
only pion and nucleon intermediate states [|19|||2^] which contribute only to the isoscalar 
polarisation. However, it was shown by Jenkins and Manohar and Jenkins ||2^ that 
it is necessary to include both the octet and decuplet of baryons in loop computations of 
the baryon axial currents, masses and weak nonleptonic decays in order to be consistent 
with experimental data. In this letter we compute the electromagnetic polarisability of the 
nucleon to leading order in chiral perturbation theory including the lowest lying decuplet of 
baryon resonances, octet of pseudo-goldstone bosons and octet of baryons as intermediate 
states. We will show that the isovector polarisability is determined at leading order by 
kaon loops, which is another example of the role the strange quark plays in the structure 
of the nucleon. 

The strong interaction dynamics of the pseudo-goldstone bosons associated with the 
breaking of SU{3) l ® SU{3)ii chiral symmetry to SU {3)v flavour symmetry are described 
to leading order in a derivative expansion by the Lagrange density 
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where /tt ~ 135 MeV is the pion decay constant and E is the exponential of the pion field 



E = exp(2zM//^ 



(2) 



with M being the octet of pseudo-goldstone bosons 

M = TT" 
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The covariant derivative is with respect to the U{1)q of electromagnetism 



D^ = d^ + tA^[Q, ] , (4) 

where A denotes the electromagnetic field. Unlike the pseudogoldstone bosons that have 
a well defined transformation under the chiral symmetry E LTiR\ the baryon fields do 
not have a unique transformation. It is convenient to define the field ^ via = E 

i = eMiM/U) , (5) 
that transforms under chiral symmetry as 

^ L^U^ = U^R^, (6) 
where U is defined implicitly by Eq. and vector fields with definite parity 

V^ = ^{^D^^^+^^D^^), A^='-{CD^^^ -^^D^^) , (7) 

where A^ transforms homogeneously and transforms inhomogeneously under chiral 
symmetry. 

As recently emphasized [|22|, processes where the baryons are nearly on their mass- 



shell, i.e. processes involving soft pions, can be consistently described by a chiral lagrangian 
involving fields of definite velocity. With such boosted fields the chiral lagrangian becomes 
a expansion in the small parameters k/ where A^ is the chiral symmetry breaking scale, 
k is the residual off-shellness of the nucleon (of order uIt^ ) and k/Ms where Mb is the 
mass of the baryon. Both of these expansion parameters are much less than unity unlike 
the expansion in P/ (of order unity) where P is the full momentum of the baryon, and 
each order in the derivative expansion is equally important. This provides a systematic 
way of power counting in 1/Mb and 1/A^ (for a review see [^). The strong interaction of 



the lowest lying octet of baryons with the pseudo-goldstone bosons is described to lowest 
order in I/Mb, lowest order in derivatives and to lowest order in the quark mass matrix 
by the lagrange density 



Cl = iTxB^{vV)B^ + 2Dl^ B^ {A^, + 2F Tr S„ S^i [A^, S„] , (8) 

where 
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and where the covariant chiral derivative is = 9^ + [V'^, ]• Only factors of the residual 
momentum are generated by derivatives acting on the baryon field. The strong interaction 
coupling constants F and D have been determined at one-loop in chiral perturbation 
theory to be -F = 0.40 ± 0.03 and D = 0.61 ± 0.04. These values correspond to an 
axial coupling of qa = 1-01 ± 0.06, much less than 1.25 extracted at tree- level. However, 
as we are computing one-loop graphs, we use the value of coupling constants extracted at 
one-loop level for consistency. 

The strong interaction of the lowest lying decuplet of baryon resonances is described 
at leading order in 1/Mb, lowest order in derivatives and to lowest order in the quark mass 
matrix by the lagrange density 



(10) 



where the elements of T„ are (we will supress the lorentz index and velocity subscript) 

2^111 _ /\++ 2^112 _ 2^122 _ A^O rp222 _ ^- j^llS _ ^ v* + 

V3 ' 73 ' ' V3 



1111 ^11) 

rpl23. _ ^ v*0 T^223 _ ^ V*- T^133 _ ^ m233 _ ^ m333 _ o- 

and Am is the mass splitting between the decuplet of resonances and the octet of baryons. 



The coupling constant Ti and C have been determined [|2| to be 7i = —1.91 ±0.7 and 
C = 1.53 ±0.3. 

The electromagnetic polarisability of the nucleon is the analogue of the Stark and 
Zeeman effects of atomic physics. For long-wavelength photons we can match onto an 
effective field theory where there are four terms at dimension seven in the chiral lagrangian 
(both isoscalar and isovector) that can contribute to the electromagnetic polarisability of 
the nucleon, 

C- D- ' ^12) 

where A, B, C, D are coefficients that are to be determined. F^'^ is the electromagnetic 
field strength tensor, N is the nucleon doublet (we only need represent the operators in 
SU{2)l ® SU{2)ii for our purposes) 

N=(l) , (13) 



and — > Ut^U^ under chiral symmetry. 

Naive power counting allows us to estimate the tree level values of the coefficients A 
through D. As they are all dimension seven operators we expect 



A ~ S ~ C ~ L> 



2A3 



O.OSGeV"^ 



(14) 



It was pointed out in that the pion cloud provides the dominant contribution to the 
polarisability and, as recently discussed in refs.|jl^[^ , the one-loop pion contribution to 
the coefficients is infrared divergent, cut off by itIt^. Therefore the one- loop contribution 
to the polarisability is enhanced by ~ nA^/rriT^ over the tree-level values. We are now in 
a position to compute the coefficients A, B, C, D to leading order in chiral perturbation 
theory from the graphs shown in fig. ^ Essentially we are integrating the pion out of the 
theory, leaving an effective theory of baryons and photons. 

Throughout our computation we will assume that the octet of baryons are degenerate, 
and the decuplet of baryon resonances are degenerate but split from the octet by Am. 
However, we explicitly retain the mass of the pion and kaon. For simplicity we write A = 
A^ + A^^ where A^ is the contribution from the octet baryons and A^^ is the contribution 
from the decuplet of resonances and similarly for B, C, D. To leading order, the coefficients 
are 

7r(D + F)2 e2 
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12 f^fUK 167r2 
which numerically are 

A' 
B' 



+ F^ + 2DF 



--D^ + F^ + 2DF 
3 



=6.5 X 10"^ GeV^ 
= - 1.4GeV-^ 
=4.4 X 10"^ GeV"^ 
= -9.7 X 10-2 GeV" 



(16) 



The isoscalar amplitude arising from the octet baryons is ~ 50 times larger than the naive 
size of the contact term. The contribution to the isovector amplitude is ~ 4 times larger 
than the naive size of the contact term. Formally this is the dominant contribution but 



6 



in reality the incalculable contact term may be significant. However, it is interesting to 
note that the isovector amplitude is dominated by kaon loops and hence strange quark 
configurations of the nucleon. 



At this point it is worth commenting on the results obtained in [|T9[. In their com- 
putation pion loops were computed in the fully relativistic theory that does not have a 
consistent power counting scheme. They found an isovector contribution proportional to 
log(m7r/M7v) which is incorrect as the nucleon mass should never appear as an argument 
of the logarithm. One might imagine infrared terms such as 1/Mb log(m7r/A^) generated 
through loop graphs at higher order in the 1/Mb expansion. However, to be consistent one 
needs to also include higher order terms in the l/A^^ expansion such as 1/A^ log(m7r/A^) 
that may be generated through loops. The only part of the computation that is unambigu- 



ous and legitimate to retain are the leading infrared -n/m.,^ terms as determined in |2C]. 
Numerically we would expect these subleading nonanalytic terms to be a few times larger 
than naive estimates of the contact terms and thus negligible. 

The contribution from the decuplet of baryon resonances is simply computed in the 
heavy baryon limit. As the formulae for the different operators are now functions of both 
the meson mass and baryon mass splitting, we introduce two functions 



R{x,y) 
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We find that the contribution to the coefficients A, B,C,D are 
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It is not sufficient to keep the leading term in a power series expansion in Am/m7r or 
Am/mx as Am, m.,^ and rriK are all comparable in magnitude and so we retain the 
complete functional form. Numerically we find that 

A^^ = 6.4 X IQ-^GeV"^ 
= -8.9 X IQ-^GeV"^ 

(19) 

C^o ^ ^ X 10-2 GeV"^ 
D^^ = 2.9 X IQ-^GeV-^ 

We see that the contribution from the decuplet adds to that from the octet for both 
isoscalar amplitudes. However, the contributions from the decuplet and octet to the isovec- 
tor amplitudes have the opposite sign and tend to cancel, but as the octet amplitude is 
much greater than that from the decuplet, this efi'ect is small. We see that it is necessary 
to include the effect of the decuplet of baryon resonances as they have a dramatic effect 
on the predictions. This should come as no surprise as their incliision is vital as pointed 
out in ref.m]. 

Recall that we treated all baryons in the octet as degenerate and all baryons in the 
decuplet as degenerate but split from the octet by Am. We have investigated the effects of 
the SU(3) breaking mass difference between the E* and the A and found that it changes 
the estimate of the isoscalar amplitudes A^^, B^^ by ~ 5% but changes the isovector 
amplitudes C^*^, D^'^ by as much as ~ 20%. However, as the isovector amplitudes induced 
by octet baryon intermediate states are much larger than those from the decuplet baryon 
intermediate states we feel confident about treating the baryon multiplets as we have. 
Other SU(3) breaking effects may be important, for instance, the coupling constants F, 
D, and C were extracted from data in the SU(3) limit. If we were to use the value of C 
from the decay A Ntv (as opposed to the average of C found from all measured decays 
T —>■ BM where T is a member of the decuplet, S is a member of the octet and M is 
a pseudogoldstone boson), we would find the contribution of the decuplet to be larger 
than our above estimates. However, to be consistent all SU(3) violation must examined at 
once, including the extraction of coupling constants F and D. Until a complete analysis of 
SU(3) breaking is performed only the couplings extracted in the SU(3) limit can be used 
consistently. 
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To compare our predictions with the experimentaUy determined values it is convenient 
to define the variables , /3^^^ and (3^^^- the isoscalar and isovector electric and 

magnetic susceptibilities. If we use the central values for F, D, C we find that 



a 



a 



(0) + ^(0) ^ 






= 13.9 X 10"^ fm^ 








= -0.16 X 10"^ fm^ 


(1)+^(1) = 






= 0.4 X 10"^ fm^ 


/3« = 




+ 


= -0.04 X 10"* fm^ 



(20) 



which are to be compared with the experimentally determined numbers 

a" = - a^^) = 12.3 ± 1.5 ± 2.0 x 10"* fm^ 

= /^(o) - = 3.5 ^ 1.5 ^ 2.0 X 10"* fm^ 
aP = + a^i) = 10.9 ± 2.2 ± 1.4 x 10"* fm^ 
(3P = /3(o) + = 3.3 ^ 2.2 t 1-4 x 10"* fm^ , 



(21) 



where the superscripts n and p denote neutron and proton respectively. These numbers 
seem to be in good agreement with our theoretical estimates , even though the experimental 
errors are still too large to be discriminating. 

In conclusion, we have computed the leading contribution to the polarisability of the 
nucleon in chiral perturbation theory. The leading infrared behaviour of the pion graphs 
had been found before [jl^[^ and contribute only to the isoscalar polarisability. To be 
consistent at the one-loop level the values for F and D were extracted from the one-loop 
computation and also the decuplet of baryon resonances was included. As expected, 



the decuplet had a large effect on the estimate from chiral perturbation theory and further, 
we have found that the leading contribution to the isovector polarisability arises from kaon 
loops and hence strange quark configurations of the nucleon. Our results are in remarkable 
agreement with the measured polarisabilities and also the model independent estimate of 
+ 13^ = 14.2 ±0.03 x 10""* fm"^ found from dispersion relations. There will be corrections 
to our results arising from terms higher order in the loop expansion, 1/Mb expansion, 
1/ expansion and from insertions of the light quark mass matrix. We expect that these 
terms are all small compared with the nonanalytic amplitudes computed in this work for 
both the isoscalar and isovector amplitudes. 



9 



The same computation can be done for the spin-dependent amphtudes. Again we 
suspect that the isovector component will be dominated by kaon loops and the baryon 
resonances will play a vital role. Similar graphs make the dominant contribution to the 
decay Ti^ A77 which proceeds predominantly through kaon loops. This will be the 
subject of future work. 

After this work was completed two preprints were released discussing the role 

of the decuplet of resonances in the polarisability of the nucleon. 
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Figure Captions 



Fig. 1. Graphs that contribute to the polarisabihty of the nucleon at leading order in 
chiral perturbation theory. We have chosen to work in v ■ e = gauge in which 
the photon does not couple directly to the baryon at leading order in 1/Mb- A 
heavy black line denotes a baryon of fixed velocity v, a thin gray line denotes a 
meson, and a wiggly line denotes a photon. 
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